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ABSTRACT 

We study the evolution of two fundamental properties of galaxy clusters: the luminosity function (LF) and 
the scaling relations between the total galaxy number N (or luminosity) and cluster mass M. Using a sample of 
27 clusters (0 < z < 0.9) with new near-IR observations and mass estimates derived from X-ray temperatures, 
in conjunction with data from the literature, we construct the largest sample for such studies to date. The 
evolution of the characteristic luminosity of the LF can be described by a passively evolving population formed 
in a single burst at z = 1.5-2. Under the assumption that the mass-temperature relation evolves self-similarly, 
and after the passive evolution is accounted for, the N-M scaling shows no signs of evolution out to z = 0.9. 
Our data provide direct constraints on halo occupation distribution models, and suggest that the way galaxies 
populate cluster-scale dark matter halos has not changed in the past 7 Gyr, in line with previous investigations. 
Subject headings: cosmology: observation - galaxies: clusters: general - galaxies: luminosity function, mass 
function - galaxies: formation - infrared: galaxies 



1. INTRODUCTION 

The subject of cluster galaxy evolution has always been 
central to cluster studies. For the early types that comprise 
the majority of the cluster galaxies, studies based on the fun- 
damental plane and/or color-magnitude relations suggest that 
their s tars are primarily old, and evolve passively with time 
(e.g. lBoweretal]ll992l IStanford et~ai]ll998t iKelson et all 

2000) . Given the hierarchical nature of structure formation, 
changes in the galaxy population over cosmic time are ex- 
pected. Numerous investigations have reported such evi- 
dence, e.g. a high er fraction of blue galaxies in more dis- 
tant c lusters (e.g. [Butcher & Oemler 1978; Ellingso n et alJ 

2001 ) , a decline in the SO galaxy abundances towards higher- 
z JPressler et all 19971) . and a population that appears to be in 
the post-starburst phase iPoggianfi et alJl999i) . 

A powerful method for quantifying galaxy populations is 
via the luminosity function (LF). The characteristic lumi- 
nosity of the cluster galaxies as a whole is reflected by 
the "knee" of the LF (M„); comparison of LFs at differ- 
ent redshifts thus will indicate the evolution of the mean 
galaxy populatio n and constrain the formation epoch of clus- 
ter galaxies (e.g. [Ba rger et all 1 19981 [De Propris et all 1 19991: 
Stra zzullo et all200fl) 

Here we present an analysis of the evolution of the cluster 
/fj-band LF from z = to z = 0.9, which has several novelties 
and advantages over previous studies. Foremost, all our clus- 
ters have mass estimates accurate to < 30% (from the X-ray 
temperature Tx). We construct the LF within the same portion 
of the cluster virial radius, which facilitates comparisons of 
clusters at different redshifts and of different masses. We also 
make use of the cluster mass information and construct the 
LF in terms of space density. Finally, our data generally are 
deeper, and cover a larger fraction of the cluster virial region. 
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These aspects allow us to integrate the LF to estimate the total 
galaxy luminosity (L) and number (A/), and study their corre- 
lation with the cluster mass (M). Based on a large sample of 
nearby clu sters, we have fo und tight correlations betw een L or 
N with M \Un et al JI2003L 12004 hereafter IIM IiM respec- 
tively), which already provide an interesting constraint on the 
regularity of the cluster galaxy formation process. 

In this study we pay special attention to the evolution of 
the N-M relation, which directly constrains simple cluster- 
ing models such as the halo o ccupation distribution (HOD; 
e.g. B erlind & Weinberg H2002) out to z ~ 1. In such models, 
the way galaxies populate dark matter halos (e.g. the number 
of cluster galaxies above some luminosity limit as a function 
of halo mass, known as the halo occupation number) are tuned 
until the observed galaxy clustering at both small and large 
scales is reproduced. Physics of galaxy formation is studied 
by identifying mechanisms that reproduce the halo occupation 
statistics (e.g. Zheng et al. 2005). Furthermore, the evolution 
of the N-M relation may revea l the mass function evolution o f 
subhalos in cluster-size halos dCoorav & Milosavli evicl20'05T) . 

In ^2] we describe the sample selection, observations, data 
reduction, and our methods for the star-galaxy separation. We 
stack the cluster data to construct the composite LFs at dif- 
ferent redshifts (50, which enable us to constrain the evolu- 
tion of and the formation epoch of the cluster galaxies. 
We then proceed to solve for the LF parameters of individual 
clusters, and examine the evolution of the L-M and N-M re- 
lations (|4}. We discuss our main findings in 50 We assume 
the cosmological parameters to be (£!m, Ho) = (0.3, 0.7, 
70/z7okms _1 Mpc -1 ). Throughout the paper we denote K s - 
band by /T-band for brevity, and adopt the Vega magnitu des 
unless noted (XAB-^Vega = 1-85, Blanton & Roweis 2006). 

2. THE DATA 

The cluster sample used here is drawn from several cat- 
alogs, and the selection is heterogeneous. All our clusters 
have measured Tx, and have been observed by either Chan- 
dra or XMM. The cluster near-IR (NIR) data were acquired 
with FLAMINGOS (hereafter FLMN) over several observing 
runs at the KPNO 4m and 2.1m telescopes, generally under 
1.1"- 1.5" seeing. Some of the basic information for the 27 
clusters in our sample is listed in Tabled including: the clus- 
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TABLE 1 
The Cluster Data 



Name 


Z 


Tx 


Klim 




Remarks" 


rt. U Cll _ \J J z, 


0355 


3 4 


19 3 


35 


2m S 


Ahell0085 


0.0551 


6. 1 


18.8 


0.39 


2m S 


Ahelll795 


0.0625 


5.5 


19.0 


0.47 


2m S 


Ahell?0?9 


0.0773 


9. 1 


19.0 


0.44 


2m S 


rt.UCllZ.Z, «_> J 


OSOA 


6 87 


18 3 


53 


2m S 


rt.UCHZ. 1 tz 


0909 


8.46 


18 7 


54 


2m S 


rt.UCllZ.Z,"'T i T 


0968 


7 1 


18 6 


62 


2m S 


Ahell2055 


0.1020 


5.8 


19.3 


0.73 


2m S 


rt.UCllZ.VJ Jt 


1 130 


7 93 


18 8 


68 


2m S 


Ahpllldl 3 

rt.UCll itl J 


1430 


6 56 


19 


89 


2m S 


rt.UCll 1 Z.VTT 


1710 


3 58 


18.9 


92 


2m S 


A hell 19 14 

.rt-UCll 1 ^ 1 *"T 


0.1712 


8.41 


18.4 


0.97 


2m S 


rt.UCll 1 U0 7 


1832 


8 58 


19 3 


1 02 


2m S 


Ahell?390 


0.228 


11.5 


19.3 


0.50 


4m J 


RYT719Q 7-+.0D0S 

ivy\. J _ 1 _ 7 . ITVUuJ 


234 


6 78 


19 4 


68 


4m J S 


Abelll835 


0.2532 


9.5 


19.3 


0.61 


4m,J,S 


Abelll758 


0.279 


6.57 


18.6 


1.70 


2mJ^S 


Abelll995 


0.319 


10.37 


19.6 


0.72 


4m,J,S 


MS 1358.4+6245 


0.328 


6.92 


20.4 


0.91 


4m,J,S 


MS 1621. 5+2640 


0.426 


7.6 


20.0 


1.09 


4m,J,S 


RXJ1347.5-1145 


0.451 


14.1 


18.8 


0.83 


4m,J 


RXJ1701. 3+6414 


0.453 


5.8 


18.7 


1.33 


4m,S 


RXJ1524.6+0957 


0.516 


5.1 


19.8 


1.59 


4m,S 


MACS0717.5+3745 


0.548 


14.4 


19.3 


0.97 


4m,J 


MACS0647.7+7015 


0.584 


11.58 


19.0 


1.16 


4m,J 


RXJ1120.1+4318 


0.6 


5.3 


20.0 


1.79 


4m,J,S 


MACS0744.8+3927 


0.686 


10.04 


19.3 


1.44 


4m,J,S 



a 2m & 4m indicate the KPNO telescope used; J means ./-band data also available; 
S means the cluster is covered by SDSS DR5. 



ter name, the redshift and Tx (in keV) of the cluster, the limit- 
ing Zf-band magnitude, the fraction of the cluster virial radius 
r2oo 7 that the FLMN field-of-view (FOV) corresponds to, and 
the available data in other bands (J and optical bands from the 
SDSS). The reduction of the data was carr ied out by a p ipeline 
developed at the University of Florida (see Elston et al. 2006). 

To increase the coverage in both redshift and mass, 
we supplement our FLMN data with those presented in 
Stanfo rd et ail J20021 hereafter IS021) . The 45 clusters in 
IS02I were not selected by any well-defined criteria, and 
the data were obtained with several different telescopes 
and instruments. These instruments have a much smaller 
FOV than FLMN. Som e of these data have been used by 
iDe Propris et alJ i 19991 hereafter Ie>P991) for the stu dy of the 
A'-band LF evolution. We include 14 clusters from S02 that 
have measured Tx and both / and /^-band data. The combined 
sample includes 41 clusters from z ~ to z ~ 0.9, all of which 
are observed out to at least r2ooo ~ 7"20o/3 (Table^. 

The varying point spread function across the FLMN array 
prevents the use of morphological information for star-galaxy 
separation from our images. For the majority of clusters at 
Z > 0.2, both J and /T-band photometry are available. In this 
case J— K > 1 can be used to select galaxies from the source 
catalog ( DP99). For the clusters that lack /-band data, we use 
data from the SDSS for help in star-galaxy separation. Ob- 
jects identified by the SDSS as point sources with iab < 21.3 

7 The virial mass is inferred from Tx using the empirical mass-temperature 
relation obtained by Finosuenov et al. 1 2001). Self-similar evolution of the 
relation is assumed. We denote the radius within which the me an overdensity 
is A times the critical density as r&. The Navarro et al. 1 1997, NFW) profile 
with concentration c = 5 is used for converting between radii of different 
overdensities (e.g. r 700 = 0.56r 20 o, looo = 0Alr 20o , r 2 ooo = 0.33r m ) We use 
a c = 3 profile for the spatial distribution of galaxies within clusters 1L04I) . 
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FIG. 1. — The evolution of apparent LFs, evaluated within r2ooo- The 
LFs are presented in terms of the physical space density (as opposed to the 
comoving density). The binning in redshift, and the number of clusters in 
each bin are: < z < 0.1 (7), 0.1 < z < 0.2 (6), 0.2 < z < 0.3 (4), 0.3 < 
z < 0.34 (8), 0.34 < z < 0.45 (4), 0.45 < z < 0.55 (3), 0.55 < z < 0.7 (4), 
0.7 < z < 0.9 (5). The midpoint redshifts are marked in each panel. 



(the 95% completeness limit for point sources) are regarded as 
stars. Based on the i—K color distribution of the sources from 
the FLAMINGOS Extragalactic Survey (lElston et alJl2006L 
hereafter FLAMEX), we find that for objects brighter than 
K= 18.06, the SDSS identifies 95% of the stars. For objects 
fainter than iab = 21 .3, we filter out the stars with the criterion 
0?-0ab > 1.4(/-/T)ab+0.99, which is a cut that we deter- 
mine will filter 90% of the stars. At K > 18.06 the galaxy 
surface density outnumbers that of the stars, and we conclude 
that our star-galaxy separation scheme will induce < 1 % con- 
tamination in the final galaxy catalogs. 

3. tf-BAND LUMINOSITY FUNCTIONS FROM z = TO z = 0.9 

The method developed in lL03l for obtaining the total galaxy 
luminosity or number of a cluster is to use the observed, 
background-corrected galaxy number and flux to solve for the 
LF parameters and assuming a fixed faint-end power- 
law slope a of the LF. Before we study the N-M and L-M 
relations at intermediate redshifts, it is important to examine 
the evolution of the LF, as this (1) allows us to understand the 
behavior of the faint-end, and (2) provides a means to apply 
the evolution correction to the observed galaxy flux. 

Using the number counts derived from the FLAMEX sur- 
vey for statistical background correction, we follow IL04I to 
construct the composite LF in the observer's frame in eight 
redshift bins (Fig. 0. For each bin, we choose a midpoint 
redshift z«,„, and adjust the ^-correction and distance modulus 
of individual clusters to best represent how the clusters would 
appear if they were all at z cm (e.g. DP99J. The combined cor- 
rections are typically less than ±0.3 mag. As we discuss in 
IL04 the main difference in the LFs between high and low 
mass clusters is the normalization 0* . Because our purpose 
here is to examine the evolution of the LF shape, we ignore 
this difference in constructing the composite LF. In each panel 
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FIG. 2. — Evolution of K t as determined within r2ooo (solid pen- 
tagons). Results from DP99 are shown as hollow circles. Other data 
points (open square s) are: two clusters at z = 0.31 jAndreon et al J I2005D . 
EIS0048 at z = 0.64 IMassarotti et all2003l) . three cluster s at z = 0.83, 0.89 
and 1.03 I Ellis & Jones 2004), two clusters at z ~ 1 I Kodama & Bower 
l200l. MG2016+1 12 at z = 1 tToft et alj|200l cluster 3C324 at z = 1.2 
1 Nakat a et all2001l) . and three clusters at z = 1.2 Jgtrazzullo et al. 20Q3). 



we show the best-fit Schechtei] (119761) function to the data, 
with a = -0.9. This value is consistent with what we find in 
L04, and has been adopted by several previous studies. We 
also note that the shape of the LFs (K* and the faint-end) does 
not change significantly within 30-60% of r2oo, which means 
the LFs shown in the Figure (evaluated within r2ooo) should be 
representative of the LF evolution. 

In Fig. 13 we show the best-fit as a function of redshift, 
and compare them with the predictions from a population syn- 
thesis model dBruzual & Charlotll200 3. hereafter BC03]). The 
model galaxies are formed with the Salpeter initial mass func- 
tion and solar metallicity in a single burst at Zf oua = 1-5, 2, 
and 3. We ha ve normalized the model predictions to the M* 
found in L04. The dot-dashed line in the Figure shows the 
with no evolution correction applied. We also s how th e re- 
sults from several previous studies, mainly that of DP99. Our 
data, together with most of the clusters at z < 1 from the liter- 
ature, seem to favor the model where the stars in clusters are 
formed at Zf orm = 1.5-2. None of the previous studies noted 
here examines the LF at a fixed fraction of the virial rad ius. I n 
the absence of any radial dependence of the M* (see e.g. lL04l) . 
studying a region that is a constant fraction of the virial region 
is not required. Because we look at the cluster galaxies as a 
whole, the behavior of K*(z) would also reflect the evolution 
of field galaxies that fall into clusters much later than the old 
spheroids. Our inferred Zf orm therefore seems low compared to 
that d erived from stud ies targeting early type cluster galaxies 
(e.g. Ellis et al.H997l) . 

4. EVOLUTION OF THE NEAR-IR SCALING RELATIONS 

The results from the previous section indicate that (1) a pas- 
sively evolving stellar population formed in a single burst at 
z = 1 .5 can describe the LF (of the whole cluster galaxy popu- 
lation) evolution based on our z < 0.9 cluster sample, and (2) 
the faint-end slope of LFs is consistent with a = -0.9 at all 
redshifts. We can then use the BC03 model prediction of the 
spe ctral aging for the k- and evolution corrections. We fol- 
low lL03l to obtain the total galaxy luminosity and number for 
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FIG. 3. — The panels on the left (right) column show the L-M (N-M) 
correlation. The top three rows are the results measured within ryoo, rioooi 
and r2ooo> respectively. Clusters in each redshift bin are plotted with different 
styles: z = 0.0 — 0.2: circle; z = 0.2 — 0.4: square; z = 0.4 -0.6: triangle; 
z = 0.6 — 0.9: star. The bottom row shows the correlations for the z < 0. 1 
local cluster sample (see text for details). The straight lines in the bottom 
two rows show the best-fit L-M and N—M correlations for the local sample. 



our clusters. Specifically, we integrate the LF over all galaxies 
more luminous thanM^>(z)+2, where Mgjz) is the restframe 
characteristic magnitude based on the BC03 model. 

We show in Fig.|3]the L-M and N-M correlations evaluated 
within 7700, Hooo, and r2ooo- Clusters are binned into four red- 
shift bins. It is interesting to see that these correlations also 
exist for the intermediate redshift clusters studied here. For 
comparison, in the bottom panels are sho wn th e correlations 
for the 93 clusters at z < 0.1 studied by IL04I (hereafter the 
local sample). For these local clusters L and N are both first 
determined within r2oo [down to Mk*(z = 0) + 21, then scaled 
to r2ooo using an NFW profile with c = 3 (e.g. LOJ)- 

We examine the existence and the form of evolution in the 
data by considering 



N(M, z)=N (l + zViM/Mof 



(1) 



for the N—M evolution, where No and Mq are the normaliza- 
tion factors of the relation. The evolution is specified by 7. 
With Mo set at 10 14 3 M Q , we fit the above expression to all 
clusters to find out the parameters (No, 7, s). It is expected 
that the parameters 7 and s would exhibit some degeneracy, 
and therefore we also include the local sample for this ex- 
ercise, whose size and the range in mass would greatly help 
determine iVo an d s. Within different regions (e.g. r7oo, nooo, 
and r2ooo)> the constraints on the evolution of the N-M rela- 
tion are consistent. For simplicity, we only present the re- 
sults from A^ooQ-^2000 correlations, which are the tightest of 
the three cases. In Fig. |4] we show the joint constraints for 7 
and s. The smaller (larger) contour is obtained when the lo- 
cal sample is included (excluded). The contours correspond 
to 68% confidence region. For each case the best-fit (7, s) 
is shown as the symbol at the center of the contour. Com- 
bining both samples results in a much tighter constraint on 
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FIG. 4. — Joint 68% confidence region for 7 and s (see Ean. HI for the 
definition). The large contour results when Eqn.[3is fit to the high-z sample. 
A tighter constraint on 7 is obtained when both high-z and local samples 
are used (the dashed contour). The points show the best-fit values for (7, 
s): (-0.53, 0.66) for the high-z sample, and (-0.03, 0.76) for the combined 
sample. 



7. If we further marginalize over s, we find that without the 
local sample, 7 = -0.53 ± 0.59; inclusion of the local sam- 
ple gives 7 = -0.03 ±0.27. Both results are consistent with a 
no-evolution scenario. 

5. DISCUSSION 

We can understand the no-evolution of the N-M relation by 
looking at the integration of the LF: 



N(M,z) = V^T(a+\Mz)/U{z)), 



(2) 



where V is the cluster volume, L; is the luminosity corre- 
sponding to the lower limit of the integration, and V is the in- 
complete gamma function. Because L; is set to be a fixed frac- 
tion of L*, it is the product of V and 4>* that determines N. For 
two clusters i and j which have the same mass but at different 
redshifts (say Zj > Zi), V) < Vj (assuming that r& oc 1 /H(z), 
where H(z) is the expansion rate of the Universe). One thus 



expects that on average, 



> 



Indeed, 0» derived for 



our clusters agrees with this expectation. That the increase 
in the galaxy number density is offset by the decrease in the 
virial volume is the main reason for the lack of evolution of 
the N-M relation. 



Previous studies have found suggestions of lack of HOD 
evolution toward z ~ 1, based on clustering st atistics 
dYan et all 120031 iPhlens et alll2006t ICojj et all 120061). In a 
recent attempt to extract both cosmological information and 
cluster physics from a large sample of clusters, it is found that 
the optical observable-mass relation, althou gh only weakly 
const rained, does not evolve with redshift ( Gladders et al. 
2006). We note, however, passive evolution has not been ac- 
counted for in this analysis (S. Majumdar 2006, private com- 
munication). Our study is complementary in the sense that 
we probe the (high mass end of) occupation number directly, 
with reliable cluster mass estimates from Tx . 

It is possible that changes in both the red and blue galaxy 
populations with redshift are hidden from the evolution of the 
N-M relation, if they conspire to make the total galaxy num- 
ber unchanged. We are in the process of acquiring color infor- 
mation in order to constrain the evolution of these populations 
separately. 

The fact that only when the large, nearby cluster sam- 
ple of L04 is included can we obtain reasonably tight con- 
straint on the N-M evolution points to the need for a 
much larger sample with reliable mass estimates at higher- 
z. This is also demanded for large optical/NIR cluster surveys 
which may rely on optical/NIR based cluster mass estimators 
(e.g. FLAMEX). Based on simple Monte Carlo simulations, 
we estimate that at least 20 clusters at z = 0.8-1 .0 are required 
to constrain 7 to 0. 1 level. While more clusters at intermedi- 
ate redshifts may provide better constraints on the evolution 
of the slope of the scaling relation, z > 1 clusters clearly will 
help pinpoint the formation epoch of the cluster galaxies. 
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